The radiative lifetime of metastable CO (a^II, = 0) 
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We present a combined experimental and theoretical study on the radiative lifetime of CO in 
the a^Ili^2,v = state. CO molecules in a beam are prepared in selected rotational levels of this 
metastable state, Stark-decelerated and electrostatically trapped. From the phosphorescence decay 
in the trap, the radiative lifetime is measured to be 2.63 ± 0.03 ms for the a^IIi, v = 0, J = 1 level. 
From spin-orbit coupling between the a^H and the A^Il state a 20% longer radiative lifetime of 
3.16 ms is calculated for this level. It is concluded that coupling to other states contributes to 
the observed phosphorescence rate of metastable CO. 



Triggered by the observation of CO in the upper at- 
mosphere of Mars with a UV spectrometer on board 
the Mariner 6 spacecraft in 1969 [1], researchers started 
studying the UV bands of CO in increasing detail. In par- 
ticular, the a^n ^ X^E+ transition between the ground 
state and the lowest electronically excited state of CO, 
the so-called Cameron bands, became the subject of theo- 
retical as well as experimental interest. From the appear- 
ance of the spectrum it was concluded that the intensity 
in this spin-forbidden transition mainly originates from 
spin-orbit mixing of the a^H state with a state; a 
detailed analysis revealed that mixing of states with 
the a^n state contributes less than 1% to the total inten- 
sity 0. Using perturbation theory and taking only the 
interaction with the A^H state into account, James cal- 
culated the radiative lifetimes for various ro- vibrational 
levels in the a^U state [3j. These lifetimes are strongly 
J-dependent, but the ratio of the lifetimes of different J- 
levels is known with spectroscopic accuracy. The J = 1 
level in the a^Hi, v = manifold, from now on indicated 
as the J) = (1,0,1) level, has the shortest lifetime, 
calculated by James as 2.93 ms; the (2,0,2) level, for in- 
stance, lives 54.66 times longer. 

An accurate experimental value for the radiative life- 
time of ro-vibrational levels of metastable CO is a bench- 
mark for theoretical calculations on spin-orbit mixing of 
electronically excited states as well as on transition dipole 
moments. It has proved to be difficult, however, to exper- 
imentally determine the lifetime of anyof the levels of the 
a^U state of CO accurately [1 d, H, These lifetimes 
can either be extracted from absorption measurements on 
the Cameron bands or, more directly, from measurements 
of the phosphorescence decay. For absorption measure- 
ments, experimental parameters like the line-integrated 
number density of ground- state CO molecules in a certain 
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ro-vibrational level, together with the absorption line- 
shape and the spectral profile of the light source have 
to be accurately known. More problematic is that the 
Franck- Condon factors for the vibrational bands of this 
transition also have to be known to be able to deduce a 
radiative lifetime. For phosphorescence decay measure- 
ments, on the other hand, the ro-vibrational state distri- 
bution of the metastable molecules needs to be known, 
and, on the timescale of the phosphorescence, collisions 
need to be avoided. To achieve this, laser excitation of 
CO to single ro-vibrational levels of the metastable state 
in the collision-free environment of a molecular beam has 
been used [8, 9, 10]. An intrinsic problem in these molec- 
ular beam experiments is, however, that the molecules 
move with a high speed, typically limiting the time dur- 
ing which the phosphorescence can be detected to a frac- 
tion of a millisecond. 

Stark deceleration of a beam of polar molecules fol- 
lowed by electrostatic trapping enables the observation 
of state-selected molecules for several seconds [11]. Re- 
cently, the radiative lifetime of vibrationally excited OH 
radicals was determined by recording the temporal de- 
cay of their population in the trap [12*]. Here we re- 
port the Stark deceleration and electrostatic trapping of 
metastable CO molecules that are laser prepared in either 
the (1,0,1) or in the (2,0,2) level. The radiative lifetimes 
are measured by monitoring the phosphorescence decay 
of the trapped molecules. Calculations to rationalize the 
observed lifetimes are presented. 

Stark deceleration of metastable CO molecules has 
been described in detail before [13]. In the present exper- 
iment, we use a Stark decelerator with 108 electric field 
stages to load the metastable CO molecules in an elec- 
trostatic trap. The experimental setup is schematically 
shown in Fig. [H A pulsed beam of CO molecules with 
a mean velocity of 320 m/s is produced by expanding 
a mixture of 20% CO in Xe from a cooled pulsed valve 
(T = 203 K). A packet of CO molecules in the upper 
(low- field seeking) A-doublet component of the (1,0,1) or 
(2,0,2) level in the metastable a^H state is prepared by di- 
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rect laser excitation from the 'u = ground state. 

After passing through a skimmer, the packet of molecules 
is slowed down in the Stark decelerator and subsequently 
loaded and confined in an electrostatic quadrupole trap. 
When in the (1,0,1) level, about 10^ CO molecules are 
trapped at a density of 10^/cm^ and at a temperature of 
around 20 mK; in the (2,0,2) level an order of magnitude 
less molecules are trapped at a somewhat higher tem- 
perature. A detailed description of the Stark decelerator 
and of the electrostatic trap that have been used can be 
found elsewhere [14] . 




FIG. 1: Scheme of the experimental setup. A pulsed beam of 
CO molecules, laser prepared in either the (1,0,1) level (exci- 
tation scheme given in the inset) or the (2,0,2) level, is slowed 
down in the Stark decelerator and stored in an electrostatic 
trap. The phosphorescence to the electronic ground state is 
measured by a photomultiplier tube (PMT). 

The ultraviolet phosphorescence back to the electronic 
ground state that escapes through a hole in one of the 
trap electrodes is imaged with a lens onto a photomulti- 
plier tube (PMT). In the upper part of Fig. [2] a semilog- 
arithmic plot of the phosphorescence signal is shown as 
a function of time. The CO molecules are prepared in 
the (1,0,1) state at t = ms, and the strong phospho- 
rescence peak around 4-5 ms results from molecules that 
pass through the trap with (more or less) the initial beam 
velocity. At 9.3 ms, the decelerated packet of molecules 
arrives in the trap center and the trap is switched on. 
After some initial oscillations, caused by the collective 
motion of the molecules in the trap, an exponentially 
decaying phosphorescence signal is observed. This sig- 
nal is shown on a linear scale in panel {b) of Fig. [21 A 
weighted least-squares single exponential fit to this data 
in the time-interval from 13 to 38 ms gives a decay time of 
2.63 ±0.03 ms, where the error is an estimate accounting 
for the systematic effect of the initial oscillations. The 
phosphorescence of the trapped molecules is measured in 
the presence of the trapping fields, i.e. in an electric field 
ranging from zero to 10 kV/cm. Under these conditions 
the opposite parity states of the A-doublet are mixed. 

For the measurement of the phosphorescence decay of 
CO molecules in the (2,0,2) level, shown in panel (c) of 
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FIG. 2: Semilogarithmic plot of the phosphorescence intensity 
from the trap region as a function of time after production 
of CO molecules in the (1,0,1) level in the source chamber 
(a). The time at which the trap is switched on is indicated 
by an arrow. Phosphorescence decay curves are shown on a 
linear scale for CO molecules in the (1,0,1) level (b) and in the 
(2,0,2) level (c). The solid lines are the results of a weighted 
least-square fit to a single exponentially decaying curve. 



Fig. [21 the repetition frequency of the molecular beam 
experiment is reduced from the normally used 10 Hz to 
2 Hz to allow for a ~ 500 ms observation time. To avoid 
detection of stray-light from the excitation laser, that 
keeps on running at 10 Hz, no phosphorescence signal is 
recorded during a 6 ms interval around the times that 
this laser is fired. A single exponential fit to the data in 
the 50 to 450 ms interval yields a decay time of 140±4 ms. 

Apart from the phosphorescence to the ground state, 
optical pumping by blackbody radiation and collisions 
with background gas can lead to a decay of the signal 
of trapped molecules; at the present densities, collisions 
between trapped molecules do not play a role. Optical 
pumping due to room temperature blackbody radiation 
is calculated to proceed at a rate of around 0.014 s~^ for 
CO in the a'^H^v = state [11]. The contribution of op- 
tical pumping to the overall trap loss rate can therefore 
safely be neglected. The loss rate due to background col- 
lisions was measured to be 0.17 s~^ for trapped ground- 
state OH and OD radicals in this apparatus. As these 
measurements were done under otherwise identical con- 
ditions, we assume the same collisional trap loss rate for 
metastable CO. For CO molecules in the (1,0,1) level 
the loss due to collisions can then also be neglected; for 
molecules in the (2,0,2) level collisions slightly contribute 
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to the observed trap loss rate. After correction for this, 
a radiative hfetime for CO molecules in the (1,0,1) level 
of 2.63 ± 0.03 ms is found, whereas the radiative life- 
time for molecules in the (2,0,2) level is determined as 
143 ± 4 ms. The ratio of these lifetimes is 54.4 ± 1.6, 
in good agreement with the ratio of 54.66 ± 0.01 that is 
expected from the ratio of the 1] = 1 character in the 
respective wavefunctions [15]. 

To calculate the radiative lifetime of the (1,0,1) level, 
we use the model as originally described by James 
Q. Essential in this model are the matrix elements 
{vA\hA,a{^)\'^a = 0) of the r-dependcut spin-orbit cou- 
pling hA.ai'f) matrix elements {vA\lJiA,x{'f)\'^x) of 
the r-dependent electronic A^Ii — X^Tj^ transition dipole 
moment /iA,x(^), together with the energy separation of 
the spin-orbit coupled states. We have computed the 
wave functions \vx)^ \^a)^ and \va)^ for the vibrational 
levels of the X^S^, the A^II, and the a^II state, respec- 
tively, with the sine- function DVR method [16]. The di- 
atomic potentials are computed with the RKR program 
of Ref. [13, taking the required spectroscopic data for 
the X-state from Ref. [HI and for the A and a states 
from Ref. [H. 

James approximated the offdiagonal spin-orbit cou- 
pling hA,a{^) as minus the diagonal spin-orbit coupling 
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In the present 



work we used the full Breit-Pauli spin-orbit operator and 
we computed the electronic wave functions at the inter- 
nally contracted multireference single and double exci- 
tation (MRCI) level employing molecular orbitals from 
a state averaged complete active space multiconfigura- 
tional self consistent field (CASSCF) calculation with 
the MOLPRO computer program. The active space con- 
sisted of the full valence space extended with a <j, a Tr^^, 
a TTy, and a 6xy orbital (in all calculations C2v point 
group symmetry was used). The Is orbitals were op- 
timized, but not correlated. The one-electron basis con- 
sists of the core- valence correlation consistent quintuple- 
zeta (cc-pCV5Z) basis set augmented with tight p func- 
tions with exponents 462.5 and 950.4 for carbon and oxy- 
gen, respectively. The exponents were energy-optimized 
in atomic calculations. The accuracy of such an approach 
was demonstrated by A. Nicklass et al. [20]. The tran- 
sition dipole moment jiiA,x{^) was also computed at the 
CASSCF+MRCI level. The same extended active space 
was used, but the orbitals for the X and A states were op- 
timized independently and the cc-pV6Z basis was used. 
The results are presented in Fig. [3l 

Using these results, we computed a radiative lifetime 
for CO molecules in the (1,0,1) state of 3.16 ms. In or- 
der to understand the 20% discrepancy with the experi- 
mental value, we performed an extensive analysis on the 
quality of the potentials, and on that of the calculations 
of the spin-orbit coupling and transition dipole moment 
functions. 

James arrived at a calculated value of 2.93 ms for the 
(1,0,1) level [3]. We repeated the calculation of James by 
employing RKR potentials based on spectroscopic data 
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FIG. 3: Panel (a) shows the A — X electronic transition dipole 
moment /J.A,x{r) (in eao) of the present work (the solid line) 
and from Ref. [21] (the dashed line), together with the diag- 
onal [ha{r)] and offdiagonal [hA,a{r)] spin-orbit coupling (in 
cm~^). Panel (b) shows the RKR-potentials of the X^E^, the 
a^^n, and the A^U states of CO and the \va — 0) vibrational 
wavef unction. 



available in 1971 and found a value of 2.94 ms. With 
our new RKR potentials and James' spin-orbit coupling 
and transition dipole function we obtain 2.93 ms. We 
therefore do not find the extreme sensitivity to the poten- 
tial reported by Sykora and Vidal, and cannot reproduce 
their value of 3.41 ms [TO]. 

We computed the spin-orbit constants of the vibra- 
tional levels in the a^II state, {va\ha{r)\va) ^ because these 
can be compared to experimental values. The computed 
values for the levels = — 8 are all between 1.4 and 1.5 
% lower than the spectroscopic values [l9|. We checked 
the dependence of the spin-orbit coupling on the one- 
electron basis set by computing hA,a and /la at r = 2.3 ao 
with a cc-pCVQZ+p and a cc-pCVTZ+p basis set. This 
gives couplings that are about 0.5%, and 1.8% higher, re- 
spectively, than the value obtained in the cc-pCV5Z+p 
basis. We determined the effect of core-valence correla- 
tion on the spin-orbit coupling in a calculation with a 
smaller active space and found that it may increase the 
couplings by about 0.5%. 

Replacing our iiA.xir) with the one computed by 
Spielfiedel et al. [2l| gives a value of 3.35 ms for the life- 
time of the (1, 0, 1) level. The ah initio method employed 
in Ref. [21] is similar to ours. The difference between 
their ijla,x{t) and ours increases at shorter distances as 
shown in Fig. [3] (a). The main reason for the difference is 
that we optimize the orbitals for the A and the X state 
separately, whereas they employed orbitals from a state 
averaged CASSCF calculation. 
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TABLE L The n^H — X^H^ transition dipole moment /in,x 
and the n^H — a^H spin-orbit coupling hn,a at r = 2.3 ao com- 
puted at the CASSCF+MRCI level (valence CAS extended 
with 27Tx and 27Ty orbitals) in an aug-cc-pVTZ basis, the en- 
ergy separation n^U — a^Il (En — Ea, with Davidson correc- 
tion) and the scaled lifetime Ts of the (1,0,1) level when the 
states 1, .., n are included. 



n 


En-Ea (cm-^) 


(eao) 


hn,a (cm~^) 


Ts (ms) 


1 


16058 


0.542 


-36.3 


3.16 


2 


48577 


-0.513 


0.1 


3.15 


3 


54284 


0.369 


-12.9 


2.74 


4 


67572 


0.016 


17.7 


2.76 


5 


70981 


0.000 


0.2 


2.76 


6 


76952 


0.303 


8.9 


2.92 


7 


81025 


-0.105 


20.8 


2.80 



We compute a total oscillator strength for the A^U^v < 
12-XiS+,v = Otransitions of 0.1764. This is 2.4% lower 
than the value of 0.1807 obtained from high resolution 
electron energy loss spectroscopy [22]. With fiA,x from 
Ref. [2l| we find a total oscillator strength of 0.1599, i.e., 
11.5% below the experimental value. 

Another test of /j.a,x is provided by the radiative life- 
time measurements of the A^H state by Field et al. [23]. 
Their deperturbed values for v = — 7 are all of the 
same quality and decrease smoothly from r{vA = 0) = 
9.9 ±0.1 ns to t{va = 7) = 8.95 ±0.1 ns. In our calcula- 
tions, however, the lifetimes increase from r(0) = 8.4 ns 
to r(7) = 8.6 ns. Other ab initio transition dipole func- 
tions, e.g., the one of Ref. [2l|, also yield lifetimes that 
increase with va- Transition dipole functions fitted to re- 
produce the measured A^H lifetimes result in computed 
lifetimes of the (1, 0, 1) level of 3.5 ms or more. 

From the results presented so far, we conclude that 
a model that takes into account only the A^H interme- 
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diate state, cannot explain the observed lifetime of the 
(1, 0, 1) level. The main sources of error are in the A — a 
spin-orbit coupling and in the A — X transition dipole 
moment. It seems unlikely, though, that the combined 
effect of these errors on the lifetime is more than 10 %. 
Only states can have both a nonzero transition dipole 
moment to the ground state and a nonzero spin-orbit cou- 
pling to the a^Hi state. Therefore, we made an estimate 
of the contributions of higher lying states by calcu- 
lating the effective a — X transition dipole moment at 
r = 2.3 ao as a function of the number of intermediate 
states, using the result to scale the value of 3.16 ms ob- 
tained above (Table [I]). The table shows that the higher 
states could account for the difference between theory 
and experiment. We suspect, however, that taking into 
account their r-dependence and obtaining convergence 
with respect to the number of states will be difficult. 

In this paper, we report the electrostatic trapping of 
metastable CO molecules. We have exploited the long 
observation time allowed by the trap to measure the ra- 
diative lifetime of two different rotational levels in the 
a^n, V = state. These two measurements are mutually 
consistent and yield an accurate value of the radiative 
lifetime of the a^Ui.v = 0, J = 1 level of 2.63 ± 0.03 ms. 
Now that the radiative lifetime of this level is known with 
unprecedented precision, discrepancies with earlier calcu- 
lated values have become apparent. This prompted us to 
perform more detailed calculations for this level. When 
only spin-orbit coupling of the a^H state with the A^H 
state is included, we compute a lifetime of 3.16 ms. We 
show that obtaining agreement between theory and ex- 
periment will require the calculation of the contribution 
of higher states. 
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